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1 Introduction 



Let d/i a (x,y) be a set of measures (in general, complex), supported on a finite set of 
products of curves in the complex x and y planes. 

Let p a , a — 2, ... ,p — 1, be a set of functions in two variables. 

Let = (x^ \..., xffl) and y^ = (y^, ■ ■ . , yffi), a = l,...,p are two sets of 
variables, where x^ and y^ are fixed by 

xf=y^=N-i, i = l,...,N (1.1) 

We shall use the following notation 

N 

» „(a+l)\ ._TT^„ 



2=1 



We consider the following integral over (2p — 2) AT variables x^ = (x^, . . . , and 
y^ = (y{ a+1 \...,y { « +1) ),a = l,...,p-l: 



a=l a=l 

where 

v 



/ n^(y (a \^)n^ a (^,^ +i )) (i.3) 

i ^ i 



^(yW.xW) = det (x^f = ]J(x^ - xf ) =: A N (x^), 



i,j=l,...,N t>3 
Jp)\ N 



Qp {y {p \ x^) = det = ]J(y^ - Vf) =■ Av(y 



i,j=l,...,N i>3 

are Vandermonde determinants, and where 



g a (y^\x^) = det(p a (yl a \xf)) , a = 2,...,p-l (1.4) 

V /i,j=l,...,N 



Developing each detp Q into iV! monomial terms (each is labeled by an element of the 
permutation group Sn), and, for given choice of the element of the permutation group, 
say a, using the change of variables inside of each iV-fold integral to the left (namely, 
X {P) _> a(x^), y(® — > cr(y^) for all (3 < a), then, using the anti-symmetry of A N (x^) 
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(which is the integrand of the very left iV-fold integral), one finds that each term of the 
mentioned development yields the same contribution. This is a standard way to re-write 
fO|) as 



1 Y 



(M) 



2p~2 



^ N {x^)\{d^\yf ) )p 2 {yf\^ ) ) 



i=l 

J i=i 

r r N 

... /nwiM.K?^) 
j j i=i 

Integrals (11.30 may be related to the so-called determinantal ensembles [6]. 

For special choice of measures dp a and functions p a , integrals ( 1 1 . 3 1> arose in the 
study of multi- matrix models, where matrices Mi, M 2 , M 3 , . . . , M 2p _2 with eigenvalues 
respectively equal to the sets {xf\ i = 1 , . . . , N} ,{yf\ i = 1, . . . , N},{x\ 2 \ i = 
1, . . . , N},. . . ,{Vi, i = 1, • • • , N}, are coupled in an open chain. It occurs in case when 
one can reduce the integration over matrix entries to the integrals over eigenvalues of 
each matrix (for these topic see [16], [17] and Appendices to [18], [1]). Depending on 
dp a and functions p a , these are models of normal matrices, and certain models of ran- 
dom Hermitian (anti-Hermitian) matrices and certain models of random unitary matrices, 
see [4], [5], [11], [3], [8], [16], [17], together with discrete versions of these matrix mod- 
els [18]. 

For instance, to obtain the partition function for the model of random N by N Her- 
mitian matrices, Mi, . . . , M 2p _2, coupled in a chain, 

2p-2 

.TrE^ 2 Vfc(M fc )+Tr( Cl M 1 M 2 +-+ C2p _3M 2p _ 3 M 2p _ 2 ) 

k=l 

one takes 

dp a {x,y) = e C2a - lxy+V2a - l{x)+V2aiy \ a = l,...,p-l (1.5) 
Pa (x J y) = e C2 « x y j a = 2,...,p-l (1.6) 
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Then x\ , i = 1,...,N, are eigenvalues of Hermitian matrices with odd numbers, say 
M 2a -i, while y\ , i = 1, . . . , N, are eigenvalues of M 2a , ot = 1, ... , p — 1. For future 
purpose, let us use the obvious freedom to re-write d\i a and p a in form 

dm(x,y) -> e cixy+yi(x) , dfip-fay) -> e"*-^- 2 ^, (1.7) 

d(i a (x,y) -> e C2Q - ia;3/ , a = 2,...,p-2, p a (x, y) -> e ^^+^)+Hi/) ) a = 2,...,p-l 

(1.8) 

In the present paper we have two tasks. 

First, we equate the integral (11.31) to the fermionic vacuum expectation value. Here 
we use the so-called p-component fermions. This may be considered as a continuation of 
of the work [11]. 

Second, as a continuation of [7], we relate Z N to the coupled p-component KP hier- 
archies, or, the same to the p component TL hierarchy. For this purpose we consider the 
following deformation of the first and the last measures 

dvi(x,y) -> d/ii(:r,y|t« n,t«) := x^e^'^+^'^d^x, y), 

dfi p ^(x,y) - d Alp _ 1 (x, 3 /|tW,n,tW) := ^e^^^" 1 ^^)^^^, y) 

OO CO 

V(x,t^) = J]i m 4 a) , V^aT 1 ,^) = ^2T m t^, a=l,p, 

m=l m=l 

and also the following deformations of functions p a , a — 2, . . . ,p — 1, 

p a (x, y) - r nW (t< Q > + [x], tW + [y]), a = 2, . . . ,p - 1, (1.12) 

where in the right hand side we have tau functions (labeled by a = 2, ... ,p — 1) of the 
one-component TL hierarchy, and where + [x] and + [y] denote the so-called Miwa shift of 
a TL (a one-component TL) higher times, details are written down below. 

The deformation f ll.9l) -( fTTT2l) relates integrals (\1.3h to the coupled p-component KP 
hierarchies. If in ( 11.31) we take the deformed measures and the deformed functions p a as 
described above, then, Zn turns out to be a certain tau function of coupled p-component 
KP, or the same, p-component TL hierarchy, where the sets of complex numbers = 
(ti , t%*\ . . . ), t*") = , . . . ), and the set of integers n^ a \ a — 1, ... ,p, play the role 

of higher p-component TL times. For the sake of brevity we shall also use the notations 
t= (t«,...,tW) andt = (t«,...,t^). 

4 



(1.9) 
(1.10) 
(1.11) 



Important to mark, that the deformation (11. 9ft . HI. 101) and f II . 121) seems do not keep the 
form ffI77l)- fO]) . In our case the interaction e c ^ M ^ M ^+i+V2 a (M 2a )+V2 a +i(M 2a+1 ) ig rep i acec i 

by arbitrary chosen one-component TL tau function (11.121) where x = is the collection 
of eigenvalues of the matrix M2 a while y = is the collection of eigenvalues of the matrix 
M 2a+1 . 

Let us note that one can consider (p — l)iV-fold integrals if he specifies the measures 
d/i a (x,y) to be proportional to Dirac delta function which equate x to a function of y (it 
may be 5(x — y)). 

The present paper is a part of series of papers devoted to fermionic approaches to 
multi-fold integrals, see [12], [1], [2]. Let us mark that our fermionic constructions of 
papers [1], [2] and of the present paper are different from what was considered in [11] and 
also different of [12]. 

1.1 Free fermions 

Let A be the complex Clifford algebra over C generated by charged free fermions {fi, fi}iez, 
satisfying the anticommutation relations 

[f h fj} + = [£, f 3 ] + = 0, [f u fj} + = 6, r (1.13) 

Any element of the linear part 

W : = (® mGZ C/ m ) © (© meZ C/ m ) (1.14) 

will be referred to as a free fermion. We also introduce the fermionic free fields 

f{x) := f(v) ■= E (1-15) 

fcez fcez 

which may be viewed as generating functions for the fj, f/s. 

This Clifford algebra has a standard Fock space representation defined as follows. 
Define the complementary, totally null (with respect to the underlying quadratic form) 
and mutually dual subspaces 

W an := (® m<0 C/ m ) © (© m > C/ m ) , W cr := (© m > C/ m ) © (© m<0 C/ m ) , (1.16) 

and consider the left and right A- modules 

F:=A/AW an , F:=W cr A\A. (1.17) 
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These are cyclic ^4-modules generated by the vectors 

|0) = 1 mod AW an , (0| = 1 mod W cr A, (1.18) 

respectively, with the properties 

/J0) = (m<0), / m |0) = (m>0), 

(0|/ m = (m>0), (0|/ m = (m<0). (1.19) 

The Focfc spaces F and F are mutually dual, with the hermitian pairing defined via the 
linear form (0||0) on A called the vacuum expectation value. This is determined by 

<0|1|0> = 1; (0|/ m / ro |0) = 1, m<0; (0|/ m / m |0) = 1, m > 0, (1.20) 
(0|/„|0) = (0|/ n |0) = (0|/ m / n |0) = (0|/ m / n |0) = 0; (0|/ m / n |0) = 0, m ? n, 

(1.21) 

together with the Wick theorem which implies, for any finite set of elements {wf. G W}, 
(0\w 1 ---w 2n+1 \0) = 0, 

(0|wi • • -W 2 n|0) = ^ S 9nv{Q\Wcr(l)Wa(2)\fy ■■■{0\Wa(2n-l)W a (2 n )\0). (1.22) 

Here a runs over permutations for which cx(l) < er(2), . . . , a(2n — 1) < er(2n) and cr(l) < 
cr(3) < ■ • • < cr(2n - 1). 

Now let {wi}i=i^N, be linear combinations of the f/s only, j G Z, and {u^}j=i,...,Ar 
linear combinations of the /j's, j G Z. Then( Ti~2"2l implies 

(0|wi- • -w N w N - ■ -wi|0) = det ((0|wiWj|0)) |ij=i,...,jv (1.23) 

Following refs. [9], [10], for all iV G Z, we also introduce the states 

(N\ := (0\C N (1.24) 

where 

C N :=h---fP_ l if AT > (1.25) 

C N :=f-i---f N if AT < (1.26) 

Ctv := 1 if iV = (1.27) 

and 

\N):=C N \0) (1.28) 
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where 

C N :=f N -x---f ifiV>0 (1.29) 

C N :=f N ---f-i ifiV<0 (1.30) 

C N :=1 if = (1.31) 

The states (11.241) and (11.281) are referred to as the left and right charged vacuum vectors, 
respectively, with charge N. 

In what follows we use the notational convention 

A N {x) = det (xf- fc )k fc =i,...,7v (N > 0), Ao(ar) = 1, A N {x) = (N < 0). (1.32) 
From the relations 

(0\f N - k f(x i )\0)=x?- k , (0|/_ JV+fc -i/(y i )iO) = yf- fe , fc = l,2,..., (1.33) 

and ( 11.23D . it follows that 

(N\f( Xl ) ■ ■ ■ f(x n )\0) = 5 n>N A N (x), NeZ, (1.34) 
(-N\f(y 1 )---f(y n )\0}=6 n , N A N (y) 1 N G Z. (1.35) 

Following [9], [10] we consider GL^ element 

g = e h , h = J2h l J l f v h i:j EC (1.36) 

Via the conjugation, (•) — > g^g^ 1 , each g G GL^ acts on the spaces (© me zC/ m ) and 
(©mezC/ m ) as linear transformations [9], [10]. 

We suppose that the following factorization condition is valid: 

9 = 9+9-, (0|2 + = (0|, <7_|0) = |0>, (1-37) 

where g+,g~ G GL^. 

Remark. Though, the property (11.3TD is valid for a rather wide class of (11. 36ft (which 
includes all cases when the sum in (I1.36P is finite) , however, we do not know the general 
theorem providing sufficient and necessary conditions to have this property in case the 
sum in (11.361) is infinite. 

Consider 

(0\v N ■ ■ -vigvi ■■■v N \0), 
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where each Vi E (© m ezC/ m ) and each Vi G (© m ezC/ m ), i = 1,...,N. Denoting Wi = 
{g + )~ 1 v i g + e (® me zCf m ) and Wi = (g^v^g^y 1 e (® me zC/ m ) we have 

(0\v N • --Vxgvx ■ ■ -v N \0) = (0\w N - ■ -WxWx ■ • -tujvlO) = det(0|u7 i w i |0)| iij= i i ... i7 v 

where the second equality is due to the Wick theorem (11.231) . Thus 

(0\v N • • -vxgvi ■ ■ - %|0) = det(0|«iflfi3j • ■ ■ |0}| iii=1) ... )JV (1.38) 

1.2 Multi-component fermions 

One obtains the so-called p-component fermion formalism by re-numerating the above 
free fermions (11.131) as follows 

fn ^ ■ fpn+a—1 j f n ^ • fpn+a—l j (1.39) 

+00 +00 

/«(*) := zk ti 0) > W*) ■= E > ( L4 °) 

k=— 00 fc=— 00 

where a = 1, . . . ,p. From (11.131) we obviously have 

[fi a \ f { I ] \ + = [fi a \ f%>]+ = 0, [fi a \ m + = 5 a ,e5 n , m . (1.41) 



Right and left vacuum vectors are respectively defined 

|(V_0> := |0), (0 i _^0| := (0| (1.42) 

p p 

where |0) and (0| were introduced in (11.190 . 
As it follows from ffTTT9l 

/W|0,...,0) = (m<0), /W|0,...,0) = (m>0), (1.43) 

(0,...,0|/W = (m>0), (0,...,0|£ Q ) = (m<0). (1.44) 

We also introduce the states 

(n( 1 ),...,nW|:=(0,0|C n(1) ---C n(p) (1.45) 

where 

C^-^fV ifnW>0 (1.46) 

C„w:= /I? •••/$, ifnW<0 (1.47) 

C n(Q) := 1 if = (1.48) 
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(1) ,..,n (p) ) :=C n(p) -- -C nW \0,0) (1.49) 



where 

Q«» :=/£,_!•••# ifnW>0 (1.50) 

C„<«o:=jS---/S ) ifnW<0 (1.51) 

:= 1 if n (a) = (1.52) 

Let us call (11.451) and ( 11.491) respectively left and right charged vacuum vectors with the 
charge (n^\ . . . , n^). 
We easily verify that 

/i a) |*,nW *> = (m<nW), ffi ) |*,nW*> = (m > «("), (1.53) 
<*,nW *|/i a) =0 (m>n^), <*,nM*|& o) = (m < n^), (1.54) 

where * serve for irrelevant components in vacuum vectors. 

Remark 1.1. For calculations we use the Wick theorem in form U.23\) . There are two 
ways to do it: 

(1) The first one is to use $1.23\) just remembering that p-component fermions are 
composed of usual ones, see M.39^l . 

(2) The second way is to use formula $1.23\) separately for each component. Namely, 
to calculate the vacuum expectation value of an operator O, first, we present it in form 

i 

Then 

(0|O|0) = 5>|Of } • • -Of |0) = $>|Of >|0> • • • (0|Oj p) |0) (1.56) 

i i 

where the Wick theorem in form $1.23\) is applied to each of (0, ... , 0\O^\0, . . . , 0). 

2 Fermionic representation for 

Consider the element of the Clifford algebra of the following form 

9 = e M g 2 e A >g 3 ■ ■ ■ e A ^g p ^e A --\ (2.1) 
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where 

Aa= I I f {a) (x)f a+1 \y)dp a (x,y), a = l,...,p-l, (2.2) 




with measure dfi a (x,y), which we do not specify. 
In (ED 

9a = e h % h a E^/^jf, ^ eC > ( 2 - 3 ) 
*.j 

so that we have 

f^9a = 9af^\ a ^(3, i6Z (2.4) 

We also suppose that each g a = e h(a \ a = 2, ... ,p — 1, may be factorized into GL^ 
elements g+ and g~ as follows (see (|1.37p ) 

a cl a a 

9a = 9a9a, (*A*\9a=(*A*\> g~ | *, 6, *) = | *, 6, *) (2.5) 

where by * we denote irrelevant components of a vacuum vector (different from the com- 
ponent a marked by hats). 

Now, let us notice that by fll .38j) we have 

(0|/ (c %i) ■ ■ • f {a) (yN)9af {a) (x N ) ■ ■ ■ f {a \x 1 )\0) = det ((01/^(^)^/^(^)10)). N 

'(2.6) 

Now let us prove, that for special choice of functions p a , a = 2, . . . ,p — 1, namely, for 

(0\f a \y)g a f^(x)\0)= Pa (y,x) (2.7) 

we have 

(My-^N, 0, . . . , 0, -N\g\0, 0, . . . , 0, 0) = Z N (2.8) 

Indeed, to get a non- vanishing expectation value in the left hand side, we have to pick 
up only iV-th term, in the Taylor series for e Al (this is because e Al is the only factor of 
g which contains the first component fermions, and the matrix element (N, *\A™\0, *) = 
until n = N). Using the known formula ( 11. 34ft . we obtain, that the left hand side of ( 12.81) 
is equal to the integral 

j dyL^J?) . . . J d» 1 (x%\yP)A N ( X W)R 1 , 
Ri = (*, 0, . . . , 0, -N\f 2 \y?) ■ ■ ■ P\yP)g 2 ■ ■ ■ |*, 0, . . . , 0) 
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where we put * on the first place of the left and right vacuum vectors to show that we 
forget about the first component fermions. This is the first step. 

A N A 

Then, we have to pick up only iV-th term, when developing the next factor e 2 . 
Otherwise, the vacuum expectation values of the second component fermions vanishes. 
This is because the second component fermions are in presence only in e Al , g 2 and in e A2 
factors of g, and the g 2 is a sum of monomials, each of which contains equal number of 
and fermions, while e Al contains only , and e A2 contains only fermions. 
Thus, second component fermions yields the expression 

f {2) (y? ] ) ■ ■■f 2) {yf)92f i2) {xf) ■ ■■f {2 \x% ) ) 

which should be integrated with measures Y\^ =1 d^(*, yf )dii{xf , *), and then substituted 
inside (N, 0, . . . , 0, -JV| and |0, 0, . . . , 0, 0). Denoting 

<0|/ (2) (yi 2) )---^ = Q,{y {2 \x^) (2.9) 

(which, by (12.61) . is equal to detp^yl , x) )) we obtain that the l.h.s of ( 12.81) is equal to 
the integral 

J d^xfj?)... j dfi 2 (x%\yP)Q3(y i3 \x^)R 2 , 

R 2 = <*,*, 0, ... , 0, -N\f®(yf) ■ ■ ■ P\yf)g 2 ■ ■ ■ |*, *, 0, . . . , 0) 

where we put * on the first and second places of the left and right vacuum vectors to show 
that we forget about the first and the second component fermions. This is the second 
step. 

Then, it is easy to see that each exponential e Aa should be replaced by their iV-th 
Taylor term, otherwise the l.h.s. of (12. 8p vanishes, it means we have 

(N,0,...,0,-N\g\0,0,...,0,0) = 

— $—i(N, 0, . . . , 0, -N\A»g 2 A% ■ ■ ■ gp-xA^O, 0, . . . , 0, 0) (2.10) 

Continuing excluding step by step third- forth- and so on component fermions, and, 
on the last step, using the known formula fjl.35j) . we obtain that (12.101) is equal to (11.31) . 

At last we want to make the following remark 
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Remark 2.1. Insert additional factors to (|2.ip as follows 

9 = e^g 2 e M g 3 ■ ■ ■ e A ^g p ^e A ^ - g := e A ^g 2 e M g z ■ ■ ■ e A *>-* g p - X e A ^ g pgi e A *> (2.11) 
where 

9a = e^ijh^f^fV, A a = J I f^(x)f a+l \y)dfi a (x,y), a = l,...,p (2.12) 

where f ip+1) (y) = f {1) (y). (Thus we add gi,g p and d/i p (x,y) to our collection of data, g a , a = 
2, . . . ,p — 1 and dfi a (x, y), a = 1, . . . ,p — 1). Then 

oo 

<o,o,...,o,o| 5o |o,o,...,o,o) = J2 CnZ n ( 2 - 13 ) 

N=0 

where cat are certain numbers and each is the following integral over 2pN variables x^ = 
(4 , . . . ,xP) and y( a ) = (y[ a \ . . . , yffl), a = l,...,p: 

Z 9= [fl e«(y (a \x (a) ) fl d^ a \y^), ^ yW (2.14) 

•* a=l a=l 

(notice that variables y^ and x^ are not fixed by (jl.lj) ) 

In (l2TTi1) dn a (x( a \ y( Q+1 )) (a = 1, ... ,p - 1) are defined by ([112]) and 

N 

dn p {x^\y^) :=l[dfi p (xf,y^), (2.15) 

i=i 

and each Q a (y^ a \ x^) is defined by (j2.6H - (j2.T[) . where now a = 1, . . . ,p. 

Sums (|2.13|) and their relation to the grand partition function of closed chains of coupled 
random matrices and to integrable equations will be considered in a forthcoming paper. 

3 Deformation of measure and relations to integrable 
hierarchies 

The described deformation 

dfii(x, y) -> d^(x, y\t (1 \ n, t (1) ) := x ni e v ^ t(a>)+v ^ x ~ 1 ' iW)) dfi 1 (x } y), (3.1) 
d/V-ifol/) - dfi p ^(x,y\&\n,t^) := y^^+^^d^-i&y) (3.2) 

oo oo 

V(x,t^) = Y,* m t\ V(x-\l^) = Y,*- m it\ a=l,p, (3.3) 

m=l m=l 
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Then, it is quite known fact that in this case Zm = tn^ 1 ', t^), where rjv(t^,t^) is 
a tau function of the (one-component) TL hierarchy. Indeed, one just re-writes (11.31) as 
2iV-fold integral with a modified measure dp mod (the latter depends on the choice of p a ) : 

r N 

^ i=l 

Moreover, as a function of t^, t^.tW, t&0, the integral Z N (t^\ &\ n u n p , t^) is 
a tau function of the coupled two-component KP, or, the same, a tau function of the 
two-component TL hierarchy, see [1]. 

Now, in addition, we consider the following deformations of functions p a , a = 2, . . . , p— 

1, 

p a {x,y)^p a {x,y\t> {a) M a \l {a) ) := (3.4) 

(n^\e H(a)(t ^f^\y[ a) ) ■ ■ ■ f a \yP)g a f {a \x^) • • • f^{x^)e R(a) ^\n^) (3.5) 

where = (t^,^, . . . ) and = (t^, ■ ■ ■ ) are the deformation parameters, and 
where the "Hamiltonians" \ k — ±1, ±2, . . . , are defined by 

oo oo +oo 

ffW(tW) = £fl< a) tW ) H^(t^) = J2 H -^k\ Hl a) =J2ti a) At ( 3 - 6 ) 

k=l k=l n=—oo 

(for future purpose we define them for a — 1, ... ,p range). 

Let us note that the expectation value (13.51) . by definition [9], [10], [15], is a tau 
function of one component TL and in our case may be denoted by 

T n(a) (t^ + [x(%iW + [ y W}) 

where 

r , fx x 2 X 3 \ 

Now let us prove that the combination of deformations (I3.ip -f l3~2l and (13.41) is equiv- 
alent to the replacement 

(JV,0,...,0,-^|0,0,...,0,0)-> 

r N (t, n, t) := (N + n« n (2) , n^, -N - n^\e H ^ge S ^\n^ , tP\ n^~ x \ -n^), 

(3.7) 



13 



where 

p oo p oo 

a=l k=l a=l k=l 

the "Hamiltonians" Hj?' were defined earlier by (13. 6p . 

Proof. Indeed, we have that each terms of type ( 12. 9p . namely, each 

(0|/ (a) (2/} Q) ) ■ ■ ■ f (a \y^)g a f {a \x[ a) ) ■ ■ • / (Q) (4 Q) )|0) = <?(v (a \s (o) ), a = 2, . . . ,p - 1, 

(3.8) 

is now replaced by 

(n ( Q )| e tf(«>(t(«>) . . . f("\yff)g a fW(x^) ■ ■ ■ fM(xff)e S(a) ^\nW) (3.9) 

which is, by definition [9], [10], [15], a tau function of one component TL. Due to ( 11.381) 
it is equal to 

detp a (xi,yj\t,n,t) 

where 

yi |t,n,t) := (nM| e ^( tCa V" (a Hyi Q) )^/ (a H4 a) ) e#(a)(i(a) V (a) > 
As for a = l,p we have 

(iV, *|e" (1) ( t(1) ) f(x 1 ) ■ ■ ■ f(x N )e Bm ^\0, *> = a 1 A JV (x)e^i ^,tW)+^ 1 ,tW) ) (3 . 10) 

(*Ae mP)(i(P)) f(y N )--J(yi)e n( ^^ (3.11) 

>.,, . k=i ft{ fc % , which contribute to the deformation respectively of d\i\ and of 
dfip-i. The end of proof. 

Thus, we obtain that the deformation of functions p a , a = 2, ... ,p — 1, and also of 
dfii, dfi p -i reduce to the fact that Z N is equal to r^(t, n, t). It is known [9], [10], [14] that 
thus constructed rjv(t,n, t) is a tau function of the coupled p-component KP hierarchy, 
or, the same, p-component TL hierarchy. 

4 Conclusion 

We equate the multi- integral (II. 3p to the fermionic expectation value (12.81) . On the one 
hand we hope that the fermionic representation allows to evaluate different magnitudes re- 
lated to the matrix models, like spectral determinants (compare with [2]), or perturbative 
series generalizing [12], [13]. On the other hand it allows to incorporate the study of these 
integrals and related multi-matrix models to the study of multi-component integrable 
hierarchies 
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